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The method of intertwining with n-dimensional (nD) linear intertwining operator £ is used 
to construct nD isospectral, stationary potentials. It has been proven that differential part of C 
is a series in Euclidean algebra generators. Integrability conditions of the consistency equations 
are investigated and the general form of a class of potentials respecting all these conditions have 
been specified for each n = 2, 3, 4, 5. The most general forms of 2D and 3-D isospectral potentials 
are considered in detail and construction of their hierarchies is exhibited. The followed approach 
" l". provides coordinate systems which make it possible to perform separation of variables and to apply 

' the known methods of supersymmetric quantum mechanics for ID systems. It has been shown that 

I in choice of coordinates and C there are a number of alternatives increasing with n that enlarge the 

^Sj ■ set of available potentials. Some salient features of higher dimensional extension as well as some 

' applications of the results are presented. 

^ ; PACS:03.65.Fd, 03.65. Ge, 02. 30.1k 

\0 \ I. INTRODUCTION 

" ^ ■ The method of intertwining provides a unified approach to constructing exactly solvable linear and nonlinear 
I problems and their hierarchies in various fields of physics and mathematics This is closely connected with the 

. supersymmetric (SUSY) methods such as Darboux's transformation, Schrodinger's factorization, and shape invariant 
'■ potential concept which deal with pairs of Hamiltonians having the same energy spectra but different eigenstates 
IO'0- general, the object of the intertwining is to construct the so called intertwining operator C which performs an 
intertwining between two given operators of the same type (differential, integral, matrix, or, operator-valued matrix 
L.J operator, etc.). In the context of quantum mechanics C is taken to be a linear differential operator which intertwines 

' two Hamiltonian operators i/o and Hi such that 
^ ■ 

PhI CHq = HiC. (1) 

. Two simple and important facts that are at the heart of the usefulness of this method can be stated as follows; (i) 
' li is an eigenfunction of Hq with eigenvalue of then ip^ = Cip'^ is an (unnormalized) eigenfunction of Hi with 
^ . the same eigenvalue Hence C transforms one solvable problem into another, (ii) When Hq and Hi are Hermitian 
O^' (on some common function space) intertwines in the other direction HqC^ = Hi and this in turn implies that 
] [Ho,C^C] = = Hi], where ^ and [,] stand for Hermitian conjugation and commutator. Therefore , two hidden 

dynamical symmetry operators of Hq and Hi are immediately constructed in terms of £ Q . These are dimension and 
form independent general properties of this method ||^. Despite this fact, like the above mentioned SUSY methods, 
the intertwining method is mostly studied in the context of one dimensional {ID) systems where £ is taken to be first 
order differential operator and Hamiltonians are in the standard potential forms. Two additional properties that arise 
in that case are that |^ ; (i) Every eigenfunction of Hq (without regard to boundary conditions or normalizability) 
can be used to generate a transformation to a new solvable problem (see Eq. (20) below), (ii) A direct connection 
to a SUSY algebra can be established by constructing a diagonal matrix Hamiltonian H = diag{HQ, Hi) and two 
nilpotent supercharges = {Q~)^ such that the only nonvanishing element of Q'^ matrix is Q21 — JC. These obey 
the defining relations of the simplest SUSY algebra 

{Q+,Q-} = H, (Q+)' = (Q-)'-0, 

which imply [H, Q^] — and emphasize in a compact algebraic form of the spectral equivalence of two ID systems. 
In the nomenclature of the SUSY quantum mechanics C is known as the supercharge operator and its zeroth-order 
(in derivatives) term as the super-potential. 
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There are important studies in the hterature which aim to generahze the SUSY methods beyond ID problems. 
These can be classified as (i) Curved-space approach (for recent studies see (l^), and (ii) Matrix-Hamiltonian 

approach [^JTl|-[r^. Both are based on the intertwining method and they mostly concentrate on extension to two 
dimensions. 

The application of the first approach to quantum mechanics was motivated by Ref. [fij which deals with free particle 



propagation on a Lie group manifold (see also |14|). Later on, this has been advanced to find the propagator of a 
free particle moving on an nD sphere Q as well as to solve both ordinary and partial differential equations with 
applications to symmetric spaces j2|. These approaches are expected to produce solvable ID n-particle problems 
from an nD free motion via some projection methods like that used in Refs. [n5ln6|. The second approach, appeared 



for the first time in Ref. [gT[, performs the extension by preserving the connection with a SUSY algebra |12| , |13| . 
This inevitably restricts the consideration to two matrix Hamiltonian such that one of them has off-diagonal entries. 
Accordingly, a matrix with elements having higher order derivative terms participates as the intertwining operator. 
This approach establishes the equivalence of two matrix systems but does not establish spectral equivalence between 
two scalar Hamiltonians. To improve it in this regard, an algorithm called the polynomial SUSY in which {Q^, Q^} 
is a polynomial of the H-matrix was introduced |13[| . 

The classification given above is by no means exhaustive; for instance one may find a nice method based on integral 
intertwining operator in Ref. (section 2.8) to generate a hierarchy of 2D problems. We should also note that 
recently the intertwining method has been used for the non-stationary Schrodinger operator [p|, pT|jT^ ]. 

The main purpose of this paper is to extend the intertwining method to an arbitrary dimension by applying it to 
a pair of nD systems characterized by Hamiltonian operators of potential form 

H, = -V^ + V,, i-0,1, (2) 

where the potentials Vi and eigenvalues of Hi are expressed in terms of 2m/'h^ and — ^^e Laplace's 

operator of i?". We shall use the Cartesian coordinates {xk\k = l,...,n}, the convention = d/dx^ and the 
abbreviation Vi = Vi{xi, throughout the paper. We purpose the ansatz that C is the most general first-order 

linear operator 

n 

£ = Lo + id = io + XI ^kdk (3) 
fc=i 

where Lq, Lk are some functions of {xk] k — 1, n} which together with Vi are to be determined from consistency 
equations of Eq. (1). In terms of the vector field L — {Li, Ln) and nD gradient operator V the operator 
Ld = ^kdk will be usually written as = • V, where "•" denotes the usual Euclidean inner product. 

In the next section by solving the first n{n+ 1)/2 consistency equations we will show that the operator is a series 
in generators of the Euclidean algebra in n-dimension. There remain n -\- 1 consistency equations which consist of n 
linear and 1 non-linear partial differential equations. Some particular solutions of these equations for an arbitrary n 
are presented in section III. In section IV we take up the integrability conditions of the remaining n linear equations 
in the context of the Frobenius integrability theory jl9j . General forms of the potentials respecting all integrability 
conditions for n = 2,3,4,5 are obtained in section V. A detailed investigation of 2D and 3-D isospectral potentials 
are given in sections VI, VII VII where we also exhibit how to generate hierarchies of potentials. 

II. INTERTWINING IN N DIMENSION: EUCLIDEAN ALGEBRA 

In view of (2) and (3) the intertwining relation (1) can be written as 

[V2 ,Ld] = - [V2 , Lo] + [Vo , Ld] + PC, (4) 
where P — Vi — Vq. At a glimpse of the right hand side of Eq. (4) and 

[y',Ld] ^J2(9^^k)dk + 2 J2idj L,)d] + 2j2idj Lk + dkL,)d,dk, (5) 

j,k j j<k 

[V2,Lo] - (V^Lo) + 2X(9,Lo)9j, (6) 

j 

[Vo, Ld] = -{LdVo) = - X L,{djVo), (7) 
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we see that the second order derivatives in Eq. (4) come, together with some first order derivatives, only from [V^, Ld\. 
Therefore by setting their coefficients to zero we obtain two sets of consistency equations: 



(9jLj=0, j = l,...,n; SjLfc + (9feLj = 0, j<k^2,...,n. (8) 

The first set gives Lj = aj + fj{x), where a^'s are constants and fj{x) depends on all of Xk's except Xj. The second 
set determines fj as fj — Cj^Xk where Cjfe's are all constants and antisymmetric in j and k : Cjk + c^j — 0. Hence 

= aj + ^CjkXk- (9) 
fc 

These solutions make the first order derivative terms at the right hand side of (5) vanish so that [V^,Ld] — 0. As 
a result of this the intertwining relation (4) simplifies to 

[y^,Lo]^[Vo,Ld]+PiLo + Ld). (10) 

From (6), (7) and (10) we get, by equating the coefficients of the first and zcroth powers of derivatives 

2d,Lo^PLf, j = l,2,...,n, (11) 
i-W^ + P)Lo^{LdVo). (12) 

These n + 1 equations constitute a reduced form of the consistency conditions for three unknown functions Lq, Vq and 
Vi. While Eq. (12) is non-linear, Eqs. (11) are linear since all components of Ld have been found. 
Eq. (12) can be considered in the following way. By virtue of 



Eqs. (11) imply that 



Combining this with (12) we arrive at 



which can be used instead of Eq. (12). 
By defining 



and using (9) Ld can be written as 



djLk^Ckj, (13) 



y'Lo^liLdP). (14) 



LoP=^Ld{Vi + Vo). (15) 



Tj = dj, Ljk = Xkdj ~ Xjdk, (16) 



Ld^Yl "J^J + ^okLjk- (17) 

3 j<k 



The generators Tj's and Ljfe's obey the following commutation relations 

K-,rfc] -0, 

[Tj,Lkm] ~ SjmTk — SjkTmi (18) 
[Ljki Lhn] = SjmLek — SjgLmk + SkeLmj — SkmLej. 

These are the defining relations of n(n + l)/2 dimensional Euclidean algebra e(n), also known as the algebra of rigid 
motion denoted by iso{n) p^ , pO| . n translational generators T^-'s form the invariant abelian subalgebra i(n) and 
n{n — l)/2 rotational generators LjkS form the semisimple subalgebra so(n). As is well known e(n) is semi-direct 
sum of t{n) and so{n) and ^ Tj = is a Casimir operator of e(n). 

Now, we shall show that the above analysis includes and naturally generalizes the well known ID case. It is evident 
that for n = 1 we have C = L{)-\- dx and P = 2L'q{x), where we take x = Xi,ai = 1 and we use the prime(s) to denote 
difFerentiation(s) with respect to the argument (when there is no risk of confusion the argument will be suppressed). 
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In that case Eq. (15) yields dx{VQ + L'q — Lq) = from which we recover the well-known forms of the ID partner 
potentials: 

Vo = Ll-L'o + b, Vi=Ll + L'o + b. (19) 

It is a standard procedure of ID SUSY quantum mechanics to take the constant b and Lq as b = Ai and Lq{x) = 
—dx\\-^4'i{x)]. When these are substituted into the first equation of (19) wc obtain : — (?!)'/(x) + Vo0i(a;) = Ai(/)i(a;), that 
is, (pi{x) is the eigenfunction of the Schrodinger's equation —(j)"{x) + Vo(j){x) = X(f>{x) corresponding to the eigenvalue 
A = Ai. Therefore the Schrodinger's equation remains covariant under the Darboux's transformations 

{cj>, Vo) ^ {L<f> = cf>'- [lncp,]'cj>, Vi=Vo- 2[\ncf>,]"). (20) 

Obviously, instead of any other fixed eigenfunction can be used to generate a transformation to another new 
potential Vi. It is this fact which allows us to apply the Darbotix's transformations successively and to construct a 
hierarchy of potentials for a given Vq- 

We conclude this section by saying that for n > 1 the differential part of the intertwining operator is a series in 
generators of e(n). In saying that we have identified the algebra generated by with e(l). A related result is that 
intertwined potentials have symmetry generators differential part of which are quadratic in the generator of e(n), that 
is, they belong to universal enveloping algebra of e(n). From now on we assume that a^-'s and Cj^s are real constants. 



III. APPLICATIONS 



Before proceeding further we consider some particular cases of Eqs. (11) and (12). 

When P = Eqs. (11) and (15) give Lq =constant and (LdVa) = 0. In view of Eq. (1) these imply, as an expected 
result, that £ is a symmetry generator of Hq — Hi : [Hq, C] ~ 0. 

Next we take P to be a constant such that P = po 7^ 0- In that case the integrability conditions djdkLo = dkdjLo 
of Eqs. (11) require that Cjk ~ for all j, k which lead to 2Lo = pod ■ r + 2b, where r = {xi, a;„) is the position 
vector and a represents the constant vector a = (ai, a„). Taking the constant b as b = a • b we get from (15) 

a - {plf+2pob-2VVo) = 0, 

which is solved by 

Vo^^py+Pob-f+gix), (21) 

where 6 is a constant vector, = J^j^] ^^'^ di^) = di^i, ■■■,Xn) is any differentiable function subjected to the 
constraint a ■ Vg{x) = 0. One may take 

ff(a;) =5(6(1) •f,...,6(„_i) -r), (22) 

such that 6(j)'s are linearly independent vectors perpendicular to a. Different choices of g define different systems 
which accept 

£t£ =_((?. V)2 + [a . (ipor + b)f - ipoa', (23) 

as a common symmetry generator. Accordingly 

££t = . v)2 + [5 . (1^0^+ b)f + ^poa^ (24) 

is a common symmetry generator for Vi = Vq +po- These also imply that Lja and £^/a are a pair of ladder operators 
for Hq: 

[i?o,£] = -PoA \Hq,C}\=pqC), \L.,C)\=pQa'. 

As a result of these we recover the existence of harmonic oscillator like spectrum in the spectrum of a class of nD 
systems described by Hq which contains many parameters and an arbitrary function. 
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Now we set all of cafe's to zero. From (11) and (15) we get Lq = f{() and P = /'(C) where / is an arbitrary 
differentiable function oi C, — a ■ r/2. Defining 

V±^^f{0±lfiO (25) 

we obtain, by virtue of (11) and (15) 

Vo = ^9{x) + V,^ ^g{x) + V+, (26) 

where g{x) may be taken as in (22). Observing that V± are form equivalent to (19) we can say that all the known 
techniques of ID SUSY quantum mechanics can equally well be used in this case. For this application the intertwining 
operator is C = f{Q + a ■ V and the symmetry generators are 

CC^ ^ a^V+ - (a ■ V)^ , C = a^V- - {a ■ . (27) 



IV. INTEGRABILITY CONDITIONS 



In this section we concentrate on the integrability conditions of n linear equations given by (11). It turns out that 
once these conditions are well understood all the consistency equations can be tackled more easily. 

By considering Lq as the (n + l)-th coordinate Xn+i = Lq of i?""*"^ and P as a function defined on it we introduce 
the 1-form 

n = dLQ- ipF, (28) 
on Here d stands for the exterior derivative and F denotes the 1-form 

n 

r = ^Ljdxj, (29) 

on i?". Now n linear equations given by (11) can be expressed as a single Pfaffian equation 17 = 0. In the Frobenius 
theory, integrability of this Pfaffian equation amounts to being able to find a positive valued integrating factor / and 
a function g such that = fdg |jl^. If this is possible then 17 = and dg = are equivalent Pfaffian equations 
and the solution (integral surface) of = is the hypersurface g ^constant. According to the Frobenius theorem a 
necessary and sufficient condition for the existence of functions g and / is the fulfillment of the so called Frobenius 
condition: 



nAdQ^O, (30) 



where A denotes the usual exterior product. 
From (28) and (29) we have 



1 " 

-[{dn+iP)dLo A F + Y,id3P)dx, A F + PdT], 



dfl 

and therefore 

n/\dn = -^[dLQ/\d{PT)-^P^T AdT], (31) 

where d in rf(PF) and dT stands for the exterior derivative of i?". The Frobenius conditions (30) is therefore equivalent 
to the following two conditions 

d{Pr) = 0, (32) 
F A rfF = 0, (33) 

provided that P ^ 0. Since both of these conditions are valid in R'\ P is defined on i?". 
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The condition (32) gives n(n — l)/2 equations 

KjkP = -2cjkP, (34) 

where 

Kjk = Ljdk - Lkdj. (35) 

Observe that Eq. (34) can also be obtained from djdkLo = dkdjLo and in deriving it we have used Eq. (13). The 
condition (33) could also be inferred from Eq. (32) upon exterior multiplication of dP A F + PdT = by F. It leads 
to n{n — l)(n — 2)/6 equations: 

LijCki] = 0, (36) 

where j < k < £ < n and the square bracket [ ] enclosing the subindexes means anti-symmetrization. Eq. (36) 
shows that any three of Lj's are linearly dependent, that is 

LjCkt + LkCij + LgCjk = 0- (37) 

Making use of (9) this can be written as 

(38) 

m 

This gives nothing in the case of n = 2 because F A dF is a 3-form and therefore identically vanishes on i?^. 

By a simple reasoning making use of the anti-symmetry of Cj^'s we see that for n = 3 the right hand side of Eq. 
(38) vanishes identically and a single condition 

L-c = a-c = (39) 

results. Here we have made use of the fact that in the case of n = 3 we have 

L = a + fxc (40) 

where c = (ci, C2, C3) = (C23, C31, C12) and "x" stands for the usual cross product of R^. For n > 4 more care is needed. 
It is not hard to check that Cjn[mCke] — for any n and hence for n = 4 the terms CiyCki], ■■■,C4^jCke] are equal to 
each other up to a sign " — " . These imply that in the case of n = 4, Eqs. (38) restrict all the coordinates to some 
constant values. But, as is evident from Eqs. (38), at the expense of constraining the form of L we can get rid of all 
these coordinate restrictions by imposing the conditions 

a[jCke] =0 ; j < k < £, (41) 
Cm[jCke]=Q ; m=l,...,n. (42) 
As is mentioned above in the case of n = 4 Eqs. (42) give only one condition 

C12C34 + C13C42 + C14C23 = 0, 

and Eqs. (41) give conditions which reduce the total number of parameters. To see this more concretely we define 
the following four vectors 

C(i) = (0, C34, -C24, C23), C(2) = (C34, 0, C41, -C31), 
C(3) = (C24, C41, 0, C12), C(4) = (C23, C31, C12, 0). 

Now Eqs. (41) can be rewritten as a - C(^j) = 0,j = 1,2,3,4. It is easy to check that, in view of Eq. (42), the 

determinant of the matrix formed by the components of the vectors C(j)'s has rank two. Therefore Eqs. (41) provide 
two of a^'s as free parameters, or, for a given a two constraints for cafe's. By taking into accoimt also (42) we get 
seven free parameters: five cafe's and two Oj's, or, three Cjfc's and four a/s. These can be chosen in many different 
ways. Moreover, one can also chose a lesser number of parameters without destroying the integrability conditions. 

For the number of conditions implied by (41-42) exceeds the number of parameters the investigation is getting 
harder and harder for n > 5. But, in the case of n = 5 one can keep again 5 of Cj^'s as free parameters by setting all 
aj's to zero. In that case Eqs. (41) disappear and Eqs. (42) give 5 constraints which reduce the number of Cj^s from 
10 to 5. Note also that, as has been done in section III, for n > 2 one can always set all Cjk's to zero and keep n Oj's 
as parameters. In such a case the condition (36) completely disappears. 

These remarks imply an important property of the intertwining method in higher dimensions; due to integrability 
conditions there are a number of choices in specifying C Evidently this fact enriches the set of intertwined potentials 
(sec the Table I in the case of n = 3). In the next section by taking cj^ =/= for at least a pair of j, k, we carry out an 
investigation which will enable us to find out the general forms of a class of potentials for n = 2, 3, 4, 5 endowed with 
mentioned richness for n > 3. 
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V. GENERAL FORM OF POTENTIALS 



By making use of Eqs. (13), (35) and (37) one can easily verify the following relations 



Z/ L 


(43) 


T 


(44) 


K C— ~) - -k^^ 


(45) 




(46) 


L ■ Vg{L^) = (2 5^ L,LjCij)g'{L^) = 0. 


(47) 



In Eq. (47) g is an arbitrary function of L"^ = L'j. Comparing Eqs. (43) and (11) we see that the general form of 
Lq is 

Lo = f{^), (48) 
provided that ^ 0. Then from any of Eqs. (11) P is found to be 

P = ^-^f'iv), (49) 

where rj — Li/Lj. Fortunately, Eqs. (44) and (45) imply that the solution (49) respects all the integrability conditions 
given by (34). 

The only equation that remained unsolved is Eq. (15) which is now as follows 

L-V{V, + Vo) = ^d,[fm. (50) 
From (46) and (47) it is evident that the general solution of this equation is of the form 

Vi + Vo = h + 2^L^, (51) 

where 2p{r])/L'^ accounts for the right hand side of (50) and h is the general solution of the homogeneous equation 

L • V/i = 0. (52) 
Hence, the general forms of Vq and Vi are, by combining (49) and (51) 



1, ^- 

-h-\ ^ 

2 L2 

1, . 



^0 = 0^+—' (53) 



^i = 2^+X^' (54) 
where 

V± = f{r,)±Ci,^f{v)- (55) 

As a result, the number of consistency equations has been reduced from (n + l)(n + 2)/2 (the sum of the munber 
of Eqs. (8), (11) and (12)) to 1, i.e., to Eq. (52). Geometrically, Eq. (52) means that at each point of the surface 
h =constant, L always lies on the local tangent space. Equivalently, L is always perpendicular to the (classical) force 
field determined by V/i. On the other hand, from group theoretical point of view Eq. (52) means that the common 
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part of the intertwined potentials is invariant under the action of the EueHdean group E{n), i.e., e^''h = h. For 
all these statements and the integrability conditions are dimension-dependent h must be determined in each case 
separately. The rest of the paper is devoted to a detailed investigation of n = 2 and n — 3 cases. 

As our investigation for an arbitrary dimension is completed two remarks are in order, (i) The above analysis 
enables us to write down a class of nD isospectral potentials provided that at least one of cafe's is different from zero. 
For instance, if only cjk then Eqs. (37) imply that L„i = for m ^ j,k and Eqs. (11) require Lq to depend 
only on Xj and Xk- In such a case, after defining tj = Lj/Lk it remains to solve Eq. (52) to find suitable n — 1 
coordinate functions, (ii) When the number of non-zero cafe's is greater than one there are a number of choices (at 
most n{n — l)/2) for rj. But, from Eq. (37) we see that these are all functionally dependent to each other. For 
example, in the case of n = 3 we have three choices t] = Li/L2,r]2 = Lx/Lz,r]z = LijLj, which obey the following 
relations 

% = 112^23 + mCsi = -Cl2. 

Instead of rji one may choose one of the variables = Li/r ■ L = Lijr ■ a, or for n = 3,ai = Li/{cxL)i. It is easy 
to verify that each of these satisfies relations similar to Eqs. (43-44) and (46) and enables us to express Lo,P,V± in 
terms of them. This freedom in the choice of coordinates once again manifests the largeness of the set of intertwined 
potentials. But, we should emphasize that these are all functionally dependent since the differential of any variable 
obeying (43) is proportional toT = L-df and therefore dT]i A dai = 0, etc. This also proves that as long as first order 
intertwining is concerned V± depend only on one variable. 



VI. 2D ISOSPECTRAL POTENTIALS 



In two dimension we have Li = (ai -|- cy) and L2 = (^2 — cx), where c = C12, a; = xi,y = X2- Prom Eq. (47) we see 
that, in terms of 

K = [L? + L2]1/2 ^ [(^^ ^ ^y-^2 ^ _ ^^)2]l/2 (5g^ 

the general solution of Eq. (52) is h = h{K), where h is an arbitrary differentiable fimction. Taking rj = L1/L2 and 
noting that L'^/L'^ = 1 + 7?^, by Eqs. (53-55) the general forms of the 2D isospectral potentials are found to be 

Vo = lh{K)+'^, V, = lh{K) + ^, (57) 

where 

V±=f{r,)±c{l+v')nv)- (58) 

In that case the intertwining operator is 

= fiv) + (ai + cy)d, + {a2 - cx)dy = f{ri) + c(l + if)d^. (59) 

As is well known, for a 2D stationary system the existence of a symmetry generator means that the system is 
completely integrable in the Liouville sense. Recalling that C and CC^ are symmetry generators of Hq and Hi, the 
potentials given by (57) are the most general forms of 2D integrable potentials which can be intertwined by a first 
order operator. 

We shall now present some examples in which for some simple forms of V- we consider the Riccati's equation (58) 
for dependent variable / and by solving it we construct the corresponding potentials. As the simplest case we take 
Vb = 0. This may happen in two different cases; (i) h = Q,V- =0, and (ii) h = — 26/k^, V- = b, where & is a constant. 
In these cases (58) is a separable equation of the form 

/2_c(l + r,2)/' = 6, (60) 

which has the general solution 

/ = {-hf/^ tan[M^(tan-i 77 - 61)] (61) 

for & < 0. This should be read as / = 6^^^ tanh[(&^/^/c)(6i - tan"^//)] for 6 > and as / = c(6i - tan^^r/)^^ for 
6 = 0, where 61 is an integration constant. Prom (58) we have 
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Vi = 2c2[k(6i - tan-i tj)]-^ , (62) 

for the case (i) and 

Vi = -2b{KC0s\- (tan"i7/-6i)]}^2. ^ -26{Kcosh[ (tan^i 7? - 6i)]}-2, (63) 

c c 

for the case (ii) corresponding to 6 < and 6 > respectively. As a result we have found a two parameter family of 
2D potentials that are intertwined to 2D free motion. Note that for 6 = —c^, 61 = we have f — crj and 

V,=2c'^ = ^^^. (64) 

As another example, taking V- = h = — and h = [2(? j k^) + 2g{K) in (57) leads us to the partner potentials 

^0 - 9{^). Vx = 5(«) + 2c2i^ (65) 

for / = cr\. In particular, for g{K) = k^, represents a 2_D isotropic displaced harmonic oscillator and Hi a 2D 
Calogero's type system for which 

Vx = 7 To + (02 - cxf + (ai + cyf. (66) 

(02 - cx)^ 

In that case for any choice of ^(k) we have £ = 0(77 + (77^ + 1)9,,]. This explicitly shows that two different families of 
potentials, such as that given by (65) can be intertwined by the same L. This is an important property that we do 
not have in one dimension. It is evident that this arises from the separability of the problem that we shall analyze 
in the next section. It is also worth mentioning that after a simple affine transformation of the coordinates and a 
restriction on (? one can easily recognize (66) as one of the four superintegrable the Smorodinsky-Winternitz 2D 
potentials [ pl| . The above particular example shows that this potential is intertwined to the harmonic oscillator and 
one of its symmetry generators is immediately obtained as LC!^ . 

VII. SEPARATION OF VARIABLES AND HIERARCHY OF 2D POTENTIALS 

The above analysis suggests the variables (k, rj) as a new coordinate system. This is a kind of the orthogonal polar 
coordinate system with displaced center in which we have 



V2 = 4{'^a,(«9«) + (1 + V^)dr,\{i + v^)dr,]}- (67) 



This implies that the eigenvalue equations of Hi accept the separation of variables in terms of {K,ri). In fact, this can 
be carried out in an easier way by introducing the coordinates 



From (59) and (67) we get 



and V2 = e-'^''P{dl + df). By defining 



/9 = — In K, ^ = — tan ^ rj. (68) 
c c 



C = .f{0 + d^, (69) 



Hp = -dl + \?^'h{p), H± = -dl + V± (e). 



and 



V± = f{0 ± /'(O- (70) 



the Hamiltonians can be written as 
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Ho = e-^^PiHp + H_), Hi = e-^^P{Hp + H+). (71) 

If we take V'°(P>0 = the eigenvalue equation Hoip^{p,^) = ^) separates as follow 

=MC/°(0, (72) 
{Hp - E°e'"'P)R{p) = -MR{p), (73) 

where M is the separation constant. For given /9-equation for Hi is the same as Eq. (73), but ^-equation is 

H_U^{S) = MU^iC). C given by (69) intertwines only solutions of to that of H+ by [/^(C) = CU^iO- 
We shall now briefly describe how to generate a hierarchy of 2D isospectral potentials. 

Taking f{i) = -<!>' {0 / (t>{0 in Eq. (70) yields V-{C} = <!>" {C} I <t>{C) ■ This is the same as Eq. (72) for M = 0. 
Therefore, each solution of (72) with M = can be used to generate a transformation to a new problem with 
potential V\. In fact, by keeping analogy with ID SUSY methods we can do more than that. For this purpose let us 
take 

y_(0=V(0-fn, /(^) = -£!§' (74) 

in Eq. (70) and suppose that the resulting stationary Schrodinger's equation 

[-dl + V{mn{^)=£nM^) (75) 

is exactly solvable, where n is a quantum number labelling the eigenvalues and eigenfunctions. If together with (74) 
we take 

/i(p) = 2e-2^qW(p)+£„], (76) 

then from Eq. (57) Vi are found to be 

Vi = 6-2^^2(^)2 + 2£^ - V(0 + n{p)]. 

In that case the separated equations of Hq are 

[-dl + V{0]U°iO = i£n + M)U^{0, (77) 

[-dl + n{p) - e-"'fE^]Rr,{p) = -{£n + M)Rn{p). (78) 

Let us choose M such that 

=£n±M (79) 

This amounts to the fact that ^-equation of Hq is the same as Eq. (75). Therefore, C/°(0 = ^n+(0 and E°,Rn{p) 
must be labelled by n+. Accordingly Eq. (78) must be rewritten as 

[-d^ + H{p) - e'-PEl^]Rn^ (p) = -£n+Rn^ {p), (80) 

The eigenvalue equation of Hi corresponding to the same E^^ can be separated such that the /9-equation is the same 
as Eq. (80) and ^-equation reads 



l-dl + 2(^)' - mWi^O = -£nUl^{0- (81) 



where 

Ui,{0 = ^U^nAO = [-|^ + 5d<^n+(0- (82) 

The function 0n(O ^hat generates the transformation is annihilated by the action of £, i.e., £(j)n{0 — 
{[?^Ti(0/?^n(C)] ~ ^£}?^n(0 = 0- Hence, in the case of M = the function ?7^(^) corresponding to 4>n{^ can 
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not be found in this way. But, by referring to a well- known theorem of the theory of ordinary differential equa- 
tions can be constructed. This theorem says that if yo{x) is a particular non-trivial solution of the equation 
ao{x)y" + a\{x)y' + a2{x)y = then the second solution yi linearly independent from yo is given by 

2/1=2/0 / (8^) 

J Vo 

Adopting this theorem to Eq. (75) where ap = — 1 and oi = the second solution linearly independent from 4>n{0 
is found to be Y{^) = 0„(O / d^/cj^liO- C generated by 0„(O applied to F(0 gives Fo(0 = ^^^(0 = -'^/KiO- 
Inserting this (as j/o) into (83) yields the desired eigenfunction corresponding to (j)n{C) 



As a result, changing the eigenfunction of Eq. (75) used to generate the transformation will lead us to a new 
eigenvalue problem given by Eq. (81). In that way a hierarchy of 2D isospectral potentials can be constructed. 

VIII. 3D ISOSPECTRAL POTENTIALS 

In order to find the general solution of Eq. (52) for n = 3 we firstly recall the intcgrability condition (39) and the 
relation (40). Secondly we observe that the set {L,c,Lxc] forms a right-handed (unnormalized) orthogonal moving 
frame which "moves" about fixed direction of c = (ci, C2, C3) = (C23, C31, 012). By using x = xi,y = X2, z = x^ we now 
introduce the variables 



/? = f ■ c, 
1^ 

Li ai + csj/ - C2Z 



1 " 1 

7= -r*. [(5+L)xc| = f-{axc) + -[{r-cf-c^r% (84) 



L2 a2 - C3X + c\z 

These obey the following relations 

(L- V)/3 = L-c = 0, (85) 
(i-V)7 = L-(Lxc) = 0, (86) 
{L-V)7i = p{ri) (87) 

where p{vi) is a quadratic polynomial in r/: 

l2 1 

Pij)) = C372 = -[(c' - + 2ciC27? + (c2 - cl)], (88) 

L2 C3 

and = cf + C2 + c^. In deriving this we assumed C3 7^ and made use of Eq. (37). 

It is now easy to see that, in view of (85) and (86), the general solution of (52) is h = h{(3, 7) where h : ^ Ris 
an arbitrary differentiable function. On the other hand, from (53-55) the general forms of the potentials are 

Vo = lh{P,j) + ^, V, = ^hil3,j) + ^, (89) 



where 



V±=f{7i)±p{7i)f{r,), (90) 



Making use of Eqs. (84-87) £ is found to be 

^ = fiv) + (5 + f xc) • V = fiv) + p{v)dr,. 
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If instead oi r] = L1/L2 had we taken one of the variables 



V2 = -r = — — — — — > m = -r 



L3 as + C2X - ciy L3 03 + C2X - ciy 

we would have obtained (L • V)r]j =Pj{r]j),j = 2,3 and V± = fjiVj) ^Pj{Vj)fj{Vj)j where 

L2 1 



^2(772) = -C2Y2 = - — [(c - 03)772 + 2ciC3r?2 + (c -Ci)], 
^3 C2 

Paim) =cij2= — [(c^ - C3)r?3 + 2c2C3r?3 + (c^ - 4)]. 

iv3 Ci 



Without any change in the /3, 7 dependence merely C would have been changed as £ = fj{r]j) + Pj{Vj)9nj- 

As an application we again consider the simplest case Vq — 0. Following an analysis similar to that made in section 
VI one can easily verify that the following 3 different potentials: 

T,(i) _ 2/f 



.^(3) = ^/.2^rL,,„,^^ 
are intertwined to 3D free motion respectively by 



A L 0\ — T] 

where /i = [&i — (cc|)~^ tan~^(p'(?7)/2c)]~^ and 61 is an integration constant. More generally, a two parameter family 
of potentials can be constructed by means of / = (—6)^/^ tan[(— 6)^/^(61 + / dri/p{'q))]. 

We shall now show that in terms of (/3, 7, if) the eigenvalue equations of Hi's accept separation of variables. Starting 
with 

dl3 = c-df, dj={Lxc)-df, dr]=^L-df, (92) 

one can easily write the differentials dx, dy, dz in terms of dj3, d'y, drj. These are as follow 

dr = + ^{Lxc)dj + ^Ldr,. (93) 

With the help of these relations the volume form dV = dx Ady Adz, the metric ds^ — dr- dr, and are found to be 

dV = -^r^-—dP A d7 A dn, (94) 
c2p(r?) 

ds' = ^2m' + ^2idl? + ^2idvf, (95) 

= c^idl + d^{L^d^)] + PMd^\p{rj)d^]. (96) 

Prom Eq. (94) we infer that the Jacobian determinant of the transformation {x,y,z) — > (/3, 7, r?) is l/(?p{rj). On 
the other hand Eq. (95) manifestly shows that the coordinate system {/3,'y,rj) is orthogonal. In virtue of (96) the 
eigenvalue equation Hq^^ {P , , i]) — ^p'^ {P , j , rj) separates, by taking tp^ {13 , j , rj) = U^{P,j)R°{r]), as 

H0,U°{(3,^)=MU°{P,j), H^R\ii) = -MR\n), (97) 

where M is a separation constant and 

= -c^L^idj + d^{L%)] + L2[l/i(/3,7) - E% (98) 
Hr, = -p{v)dMv)dr,] + fiv) - P{v)f'{v)- (99) 
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At this point we will be content with saying that by following the similar steps as for section VII one can construct 
hierarchy of 3-D isospectral potentials. 

Finally we would like to emphasize that the 3D potentials we have found depend on six parameters such that a 
large number of potentials can be generated by setting some of them to zero, or, to some particular values. Possible 
choices of parameters are represented in the Table I. The corresponding potentials can be read off from the expressions 
in the main text. 



IX. CONCLUDING REMARKS 



Main results of this study can be summarized as follows. We have studied a pair of nD Hamiltonians of potential 
forms that intertwine by first order operator C and proved that the differential part of C is an element of the 
Euclidean algebra e(n). These imply that so-intertwined systems have symmetry operators whose differential parts 
belong to enveloping algebra of e(n). The integrability conditions of consistency equations are dimension dependent 
and therefore have been considered for each case separately. The general form of potentials have been specified for 
n = 2,3,4,5 where only one linear partial differential equation which determines the common part of the potential 
remains unsolved. We have found the general solution of this equation in cases of n = 2 and n — 3. 

Three distinctive features of the higher dimensional extension of the intertwining method are that: (i) The method 
suggests coordinate systems which allows us to do the separation of variables and to utilize, in one of the variable, all 
the methods of the ID SUSY quantum mechanics, (ii) In the choice of this variable and C itself one has a number 
of alternatives increasing with n. This fact enlarges the set of available potentials for each n > 3. (iii) There exist 
families of potentials accepting the same intertwining operator. 

2D and 31? isospectral potentials we have obtained involve two arbitrary functions. The former constitute the 
most general integrable potentials which admit first order intertwining. Particular forms of these potentials may be of 
special interest for various purposes. Having in mind the projection techniques which produce exactly solvable lower 
dimensional problems from the higher dimensional one-particle problems |l5| , [l^ our analysis in section VI- VIII must 
be continued also for n — A and n = 5. As is implied by the last example of section VI, it seems to be possible to 
investigate connections among the superintegrable potentials |^ as well as to construct related potentials by repeated 
Darboux's transformations in the context of intertwining method. Velocity dependent, stationary and non-stationary 
problems can as well be considered within our approach. Work on 2D and 3D isospectral potentials which are 

at the same time superintegrable is in progress. 
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TABLE L In 3-dimension the special choices of parameters and corresponding coordinates. Note that in each case further 
choices are possible. For example, in the first three cases Cj which does not appear in the first column can be set to zero and 
instead of ri, one can also use 772, or, 773. As a completely different case in which all Cj^'s are zero has been presented in section 
III for any n > 1. 







27 = - 


V 


pi-n) 


ai — 0; a2C2 + a^c^ = 
a2 — 0; flici + a3C3 = 
as = 0; aici + a2C2 = 


f- c 
r • c 
r • c 


2f ■ (axe) + if - c)'^ - c^r^ 
2f- {axe} + if - c)^ - c^r"^ 
2r ■ (axe) + (f • cf - c^r^ 


a2 — C3X + ciz 
il+C3a-C2Z 

— C3X + C1Z 
ai+C3j;-C2Z 
a9 — C3 3; + ci z 


p{n) 
p{v) 
piv) 


ai = = 02; C3 = 
ai = = a3; C2 = 
a2 = = 03; ci = 


ClX + C2y 

ClX + C32 

C2y + Ciz 


2a3(eiy - e2a;) + {cix + C2yY - {cf + ci)r^ 
2a2{c3X - ciz) + (cix + c-j,z)'^ - (e? + ci)r^ 
2ai(c2« - C3y) + (e2y + C3z)^ - {c\ + ci)r^ 


— C22 

03 + 021 — ciH 

C3«/ 
a2 — C3X + c-iz 
ai+C3y-C2Z 

— Cr^X 


P2{V2) 

p{v) 
p{v) 


tsi = 0; C2 = = C3 
a2 = 0; ci = — C3 
03 = 0; ci = = C2 


ClX 

C2y 
C3Z 


2ci(o22/ - a3z) - ci{y'' + z'') 
2c2{—a3X + aiz) — C2{x'^ + z'^) 
203(022; - aiy) - cl{x'^ + y^) 


a2 + ciz 

ai—C2Z 
a^ + C2X 
a-i+c^y 
a9 — C3a; 


P3{V3) 
^2(772) 

p{v) 


ai = a2 = 03 = 


r ■ c 


[r ■ c) ~ c r 


C31/-C2Z 

— c^x + ciz 


p{v) 
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